DYNAMICS OF THE THERMOHALINE CIRCULATION 
UNDER UNCERTAINTY 



WEI WANG, JIANHUA SUN, AND JINQIAO DUAN 

Abstract. The ocean thermohaline circulation under uncertainty is 
investigated by a random dynamical systems approach. It is shown that 
the asymptotic dynamics of the thermohaline circulation is described by 
a random attractor and by a system with finite degrees of freedom. 



1. Introduction 

The ocean thermohaline circulation(THC) involves water masses sinking at 
high latitudes and upwelling at lower latitudes. The process is maintained 
by water density contrasts in the ocean, which themselves are created by 
atmospheric forcing, namely, heat and freshwater exchange via evaporation 
and precipitation at the air-sea interface. Thus the ocean THC is driven 
by fluxes of heat and freshwater through the air-sea interface. During the 
THC, water masses carry heat or cold around the globe. Thus, it is believed 
that the global ocean THC plays an important role in the climate |24j . 
The formulation and analysis of mathematical models is central to the 
progress of better understanding of the THC dynamics and its impact on 
climate change. Apart from a detailed modeling of the climate system us- 
ing coupled general circulation models, sometimes simplified climate mod- 
els could give insight into the general characteristics of the climate system. 
The most simplified climate models can be described in terms of stochastic 
differential equations [H)J IH HZj- These stochastic climate models can be 
viewed as comprehensive paradigms or metaphors for particular features of 
the climate system. 

We consider a two-dimensional thermohaline ocean circulation model in the 
latitude-depth (meridional) plane, in terms of the stochastic Navier-Stokes 
fluid equations (vorticity form) and the transport equations for heat and 
salinity, together with air-sea flux or Neumann boundary conditions. The 
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noise in the Navier-Stokes equations is due to various fluctuations such as 
random wind stress forcing. 

We intend to investigate the characteristics of the THC's dynamics when 
some random effect is taken into account. Our approach here is a random 
dynamical systems approach 

This paper is organized as follows. In the next section we present the THC 
model, and discuss the well-posedness of this model in Section 3. Section 
4 is devoted to the investigation of the dynamical behavior of this model: 
random attractor and finite dimensionality. 



2. A Model for the Thermohaline Circulation 

We consider the ocean thermohaline circulation in a bounded domain, i.e., 
a square 

D = {(y, z) : -/ < y < I, < z < d}, 

on the meridional, latitude-depth (y, z) -plane, as used by various authors 
j2Hl EQ I2H 03 U2j • It is composed of the Boussinesq version of the Navier- 
Stokes equations for oceanic fluid velocity (v(y, z,t),w(y, z,t)) and trans- 
port equations for the oceanic salinity S(y, z, t) and the oceanic temperature 
T(y, z, t) in dimensional form: 

v t + vv y + wv z = —p y + uAv + noise , 

w t + vw y + ww z = —p z — g(asS — axT) + i/Aw + noise , 

v y + w z = 0, 

T t + vT y + wT z = k t AT, 

S t + vS y + wS z = k s AS, 



where as and are the coefficients of volume expansion for salt and heat, 
respectively; g is the gravitational acceleration; v is the viscosity; and k$ 
and kt are salt and heat diffusivities, respectively. The density is p = 
Po(l + ozsS — axT) with p the mean sea water density. The noise in the 
Navier-Stokes equations is due to various fluctuations such as random wind 
stress forcing [THl 123 UH] • Presumably, the noise also affects the transport 
of heat and salinity to some extent, but we will ignore this effect. 
As discussed in j2Z| |2U|, this may be regarded as a zonally averaged model 
of the world ocean. The effect of the rotation can be parameterized in the 
magnitude of the viscosity and diffusivity terms. Introducing the stream 
function i/i(y, z, t) for the velocity field, 



v = -ip z , w = ip, 
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we can rewrite the above model in the vorticity form with only three un- 
knowns ip,T, S: 

(1) Aip t + Jty, Aiff) = g{a T T y - a s S y ) + z/A 2 ^ + W x , 

(2) T t + J(^,T) = k t AT } 

(3) S t + J(tP,S) = n s AS, 

where Wi(y,z,t) is a Wiener process defined on a underlying probability 
space (fl, J 7 , P) to be specified below. The noise is described by the gener- 
alized time-derivative of the Wiener process. The fluctuating noise in the 
oceanic fluid equation is usually of a shorter time scale than the response 
time scale of the large scale oceanic THC We thus assume the noise is white 
in time (uncorrelated in time) but it is allowed to be colored in space, i.e., 
it may be correlated in space variables. Note that the pressure field p is 
eliminated in the vorticity form of the Navier-Stokes equations. 
Boundary conditions for the oceanic fluid are no normal flow and free-slip 
(no-stress) on the whole boundary [19J: 

(4) ip = Atj) = on dD. 

The boundary conditions for temperature and salinity are Neumann type. 
On the air-sea interface z = d, the heat/temperature flux and freshawa- 
ter /salinity flux are prescribed as 

(5) T z (y,d) = \(6(y)-T), S z (y,d) = F(y). 

Here 8(y) is the prescribed (known) atmosphere surface temperature and 
F(y) is the mean freshwater flux (known). Moreover, A = po ^ T KT , with Br 
being the surface exchange coefficient of heat and C p the heat capacity. 
Zero flux boundary conditions are taken for T and S on fluid bottom (z = 0) 
and on fluid side (y = ±/): 

(6) T t {y,0) = 0, S z (y,0) = 0, 

(7) T y (±l,z) = 0, S y {±l, z) = 0. 

The THC model above involves stochastic and deterministic partial differ- 
ential equations (PDEs) and Neumann boundary conditions. 

3. Well-Posedness 

In this section we will show that (l)-(3) defines a well-posed model. First 
we introduce some function spaces from the theory of partial differential 
equations. 

Let W^iD) be the Sobolev space of function on D with the first generalized 
derivative in L2(D), the function space of square integrable functions on D 
with norm and inner product 

\\ u \\l 2 = ( \u(x)\ 2 dD)2 , (u,v)l 2 = / u(x)v(x)dD, u,v G L 2 (D) 

J D J D 
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The space W 2 {D) is equipped with the norm 

\\ u \\w% = \ \ u \\l 2 + \ \d y u\\ L2 + ||<9 z m||l 2 

o 

Let the W\{D) be the space of functions vanishing on the boundary dD of 
D. The norm of this space is defined as 

( 8 ) IMIw^(d) = \\9v u \\l 2 + \\d z u\\ L2 - 

Similarly, we can define function spaces on the interval (0, d) denoted by 
L 2 (0,d) and W£(0,d). 

Another Sobolev space is W^{D) which is a subspace of W 2 l (D) consisting 
of functions u with zero mean: J D udD = 0. A norm equivalent to the 
H^-norm on W 2 l (D) is given by the right hand side of (8). For functions in 
L 2 (D) having the same property, we write as L 2 (D). 
There exists a continuous trace operator. 

ldD :W 2 \D)^m{dD). 

Here H^^dD) is a boundary space, see Adams [l j Chapter 7 or below. 
Similarly, we can introduce trace operators that map onto a part of the 
boundary of 3D for instance for the subset {(y, z) e D\z = d} of D. For 
this mapping we will write 

(9) lz=d :W^(D)^H^0,d). 
The adjoint operator 

1 : =d :(H 1 H0 } d))'^(W 2 1 (D)y 

is also continuous, where ' denotes the dual space for a given Banach space. 
If set vorticity q = Aip we can homogenize boundary conditions to obtain: 

(10) q t + J(7p,q) = g(a T T y - a s S y ) + vAq + Wi , 

(11) T t + J(^,T) = KT AT + 7 ; =d (A(%)- 7z = rf T)), 

(12) S t + J(^,S) = KsAS + jUjFiy). 
New homogeneous boundary conditions are: 

ip = 0, q = on 3D, 
T z (y,d) = 0, S z (y,d) = 0, 
T 2 (y,0) = 0, S z (y,0) = 0, 

Ty(±l,Z)=0, Sy(±l,Z)=0. 

For convenience, we introduce the vector notation for unknown geophysical 
quantities 

(13) u = (q,T,S). 

Now we can define the linear differential operator from (10)-(12) 
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-uAq 
Au= | -k t AT 

We assume that F and 9 G L 2 (0, d). Note that 
d f f d 

— J Sdydz = J F(y)dy = constant. 
It is reasonable (see [T2]) to assume that 

f F(y) = 
Jo 

and thus J D Sdydz is constant in time and we may assume it is zero: 

Sdydz = 0. 

D 

Thus we have the usual Poincare inequality for S. However, this is not true 
for T. Fortunately we can derive the following Poincare inequality for T 

(14) ||T|| 2 < C (Q)(|| 72=1 T||i 2 + ||VT||i 2 ), 

as in Temam [2H] where c(O) is a constant dependent on Q. 

Introduce the phase space for our system H = L 2 (D) x L 2 (D) x L 2 (D) with 

the usual L 2 inner product and V = W\{D) x W 2 (D) x W 2 (D). 

It is obvious that the linear operator A : V — > V is positive definite. And 

define the nonlinear operator G(u) := G\{u) + G 2 (u) where 

/ -J(M 
G 1 (u)[y,z]= -jy>,T) I [y,z]- 
\ -J(i>,S) 

and 

(g{a T T y - a s S y ) 
lUimv)-l z =dT)) | [y,z\. 
iunv) 

Then the THC system can be rewritten as a stochastic differential equation 
on V: 

du 

(15) — + Au = G(u) + W u(0) =u eH, 
dt 

where W = (Wi,0, 0), W is a white noise as the generalized temporal 
derivative of a Wiener process W with continuous trajectories on M and 
with values in L 2 (D). It is sufficient for this regularity that the trace of 
the covariance is finite with respect to the space L 2 (D): tr^ 2 Q < oo. In 
particular, we can choose the canonical probability space (O, J 7 , P) where 
the set of elementary events Q consists of the paths of W and the probability 
measure P is the Wiener measure with respect to covariance Q. 
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Through integration by parts and direct estimation or from [0] we have the 
following lemmas. 

Lemma 3.1. The operator G\ : V — » H is continuous. In particular, we 
have 

(G 1 (u),u) = 0. 

and 

Lemma 3.2. The following estimation holds 

\\G 2 (u)\\ v > < ci|H|y + c 2 . 
for some positive constants c\, c 2 . 

In the following we need a stationary Ornstein-Uhlenbeck process solving 
the linear stochastic equation on D 

(16) ^ - v{k + 1)A V = m 

with the homogeneous Neumann boundary condition at 3D. Here k can be 
seen as a very large controlling parameter. 

Lemma 3.3. Suppose that the covariance Q has a finite trace : tr£ 2 Q < oo. 
Then (16) has a unique stationary solution generated by 

(t,uj) -> rj(9 t u). 

Moreover, t] is a tempered random variable in W^iD) and has trajectories 
in the space L z 2 oc (R; W^D)). Then Z(u) = (r)(uj), 0, 0) is a random variable 
in V . 

For the proof we refer to jH]. 
If we set 

(17) (q, T,S)=v:=u-Z=(q-r ] ,T,S), 
then we obtain a random differential equation in V 

(18) -n + Av = G 1 (v) +G 2 (v + Z(9 t u)), v(0) = v e H, 
at 

where G 2 = G 2 + (J(rj, q) + J(A _1 g, Arj) + J(r), Ar/) - ukArj, 0, 0). 

The above equation f)18|) is a differential equation with random coefficients 
then it can be treated sample- wise for any sample uj. We are looking for 
solution v in 

C([0,r];H)nL 2 (0,r;V), 

for all r > 0. If we can solve this equation then u := v + Z defines a 
solution version of (15). For the well posedness of the problem we now have 
the following result. 



STOCHASTIC THERMOHALINE CIRCULATION 7 

Theorem 3.4. (Well-Posedness^) For any time t > 0, there exists a 
unique solution of n~8\) in C([0, t]',H)D L 2 (0, t; V). In particular, the solu- 
tion mapping 

R + x n x H 3 [t, w, v ) -> v(t) G H 

is measurable in its arguments and the solution mapping H 3 Vq — > v(t) G H 
is continuous. 

Proof. By the properties of A and G\ (see Lemma 3.1), the random differ- 
ential equation (18) is essentially similar to the 2-dimensional Navier Stokes 
equation. Note that G2 is only an affine mapping. Hence we have existence 
and uniqueness and the above regularity assertions. □ 

Now we can define a random dynamical system since the solution mapping 

R + x Q x H 3 (t,u,v ) — > v(t,uj,v ) =: tp(t,(J,v Q ) G H. 

is well defined. First we define a so-called metric dynamical system (fi, J 7 , P, 
(6 t )teR). t : Vt — > Q, t G R is a family of measure preserving transformations 
such that (t, a>) 1— > ^u; is measurable, 9 = id,9 t + s = 9 t 9 s , for all s,t G R. 
Furthermore, the shift t is ergodic if we define it as 

w(-, 9 t uj) = w(- + t, uS) — w(t, to) for t G R 

which is called the Wiener shift. A random dynamical system (p(t,u,u) is 
well defined now since the cocycle property 

<f(t + r, u, u) = (p(t, 9 T u, <f(r, lu, u)) for t, r > 
y?(0, u,u) = u 

for any w G fl and u E H . For more detail about random dynamical sys- 
tems we refer to [2]. 



4. Random dynamics 

In this section, we investigate random dynamics of the THC First we are 
going to show that the random THC model is dissipative, i.e., it has an 
random absorbing set in the following sense: 

Definition 4.1. A random set B = {B^uj)}^^ consisting of closed bounded 
sets B{uj) is called absorbing for a random dynamical system ip if we have 
for any random set D = {Diu)}^^, D{oj) G H bounded, such that t — > 
sn P y eD(8 t Lu) \\v\\h has a subexponential growth for t — > ±00 

<p(t,u,D(u)) C B{9 t uj) fort>t (D,u) 
( l<)) <p(t,0- t u},D{e-tU})) C B(uj) fort>t (D,u). 

B is called forward invariant if 

<f(t, u, Uq) G B(9 t u)) ifuoEB{ui) fort>0. 
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Consider THC system separately. For this we introduce the spaces 

H = L 2 {D) x L 2 {D), 
V = W 2 \D) x W 2 \D). 

We choose a subset of dynamical variables of our system (10)-(12) 

(20) v = (T,S). 

Applying the chain rule to \\v\\ 2 H , we obtain by Lemma 3.1 

d 



(21) di M " + 2 ^||VT||i 2 + 2 K5 ||V5||i 2 

=2\(6(y), lz=d T) L2 - 2X( lz=d T, lz=d T) L2 + 2(F(y), lx=d S) La . 
Cauchy-Schwarz inequality yields the following estimates 

2X(9(y) nz=d T) L2 -2\{ lz=d T, lz=d T) L2 < -||0(y)||| 2 + (a - 2)\\\ lz=d T\\l 2 , 

(X 

where a is a positive constant that satisfies 2 > a > 2 — Applying the 
trace theorem || -f z=d S \\\ 2 < c 3 1| -y^ =rf S'|| ^ i < c 3 ||S , || w -i and for any e > we 
can find a c 4 (e) such that 



2(F(y), lz=d S) L2 <e\\S\\^+c,(e)\\F(y)\\l 2 , 



for some e. 

Thus by using the Poincare inequality for S G W 2 (D), (14) and choosing e 
small enough, we conclude 

(22) j t \Ml + ^\\vn\ + \M\)<^ 

where a, c$ is determined by A, kt, ks, \\@\\l 2 , 1 1 F(y) an d is defined 
by (V,, : V.V ?/ .,S). 

By the Gronwall inequality, we finally conclude that 

(23) \M% < \\m\\%e- at + °X 

Then we can easily see that the ball B(0, i? x ), where R\ = absorbs v in 
the sense of definition 4.1. 

To prove the dissipativity of the dynamical system ip we have to study ||?||z, 2 . 
From (18), we have 

(24) q t = -J(ip, q) + g(a T T y - a s S y ) + uAq - ukArj. 
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Then 

(25) —HtflL =(-m q),v)~ HI Vg||L - vk(Ar), q) 

+ g(a T T y - a s S y ,q). 

From the definition of «/(•••) and the Cauchy-Schwarz inequality we have 

(J(^,q),v) = (J(M,q) < l|Vg||L 2 ||g|kl|V77|| La 
<^Q\\1 + ^V\\IMI 

IS . | . | n ^|| i|9 ^ II ||9 - - . | - *~> I, ■ | a 

< ^h\\wi + + ^IMWMU + 

o 

Here c 6 is the constant in the Poincare inequality for rj EW^ (0, d). And for 
any e > we can find c 7 (e) > and c 8 (e) > such that 



0(a T Tj, - a s S y ,q) < #11 Vv|y|g|| L2 

4 



<j\\q\\l 2 + C7(e)~g 2 \\Vv\\% 



—vk(Arj,q) = vk(Vr),Vq) 

where g is determined by oct, ots and g, Ai is the first eigenvalue of the 
operator —A on (0, d) with Neumann boundary condition. 
Collecting all these estimates, we have 

<l lg||l 2 <- 7 (MlklL+^)||V^||| + r(M 



dt 

< - 7 (MlklL + r(M + - ^^\\v 



, M j , , , S(e)c 5 _ 5(e) d 2 

a a dt H ' 



where 7(0;) = Aiz/ - e - £|M|^i, 5(e) = c 7 (e)g 2 and 



(26) r(u) = (2\ lC& (e)k 2 v 2 + u)\\ V \\ 2 w i + giMfe- 



10 WEI WANG, JIANHUA SUN, AND JINQIAO DUAN 



Then Gronwall inequality yields 
(27) 

||?(t,w,« )||L<||«b||| r e-J'o7(«.«)* 



a a as a 



< \\v \\ 2 ™p~fo'r( e * u ) da 



r-t 

+ ' 







•^^ + r(M)e- / * 7(SrW) *"da 
1 a 

+ 6 M f I|t5||2 7 (0 aa ,) e -/.S«»r U )«tr da + ^£1^(0)11^ i e -/ S(^. 

We will show that the right hand of (27) is finite as t — > oo. In fact we have 
Lemma 4.2. // £/ie controlling parameter k is large enough that 

Ai > 



tr L , 2 Q 



and e < ^ then 



E-f(uj) > 0. 

2 



Proof, ltd formula applied to yields 

1117(^)11^+2(^+1)1/ / 11^(^)11^^=11^(^)11^+2 / fadtiO^+rtr^Q. 



Hence we can easily get that S||j7||^-i < ^h- Then Ej(u) > 0. □ 
Now we can estimate the ||?||x 2 - First we have 

limdKHlr + —(c 5 t + \\v(0)\\%))e- f ^ {6sU))ds = 0, Ra.s. 

And note that ||t>||^ is bounded by ||v(0)||^.e -£rf + ^ which tends to ^ 
exponentially We replace u by 6L t co> to construct the radius of the absorbing 
set. Then we have 

hm [ (r(9^ s u) + ^||^|| 2 - 7 (^_ s a;))e-/^^) d ^5 



= lim F(r(0su) + ^||5||i7(M)e-^ 7( ^ )dr cis 

=:i2|(a;) < °°- 
Collecting all the estimates we have 

Lemma 4.3. The random set {5(u;)} we n given by closed balls B(0, R(u)) 
in H with center zero and radius R 2 (uj) := R\ + R\{oj) is an absorbing and 
forward invariant set for the random dynamical system p generated by Alty) . 
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For the application in the following we need the particular regularity of the 
absorbing set. To this end we introduce the function space 

H s :={ueH : \\u\\ 2 s := \\A?u\\ 2 H < 00} 

where s G ffi. The operator A s is the s-th power of the positive and sym- 
metric operator A. Note that these spaces are embedded into the Sobolev 
spaces H a , s > 0. The norm of these spaces is denoted by || • \\h»- 

For our aim we show that v(l,u,D) is a bounded set in 7i s for some s > 0. 
Consider t||f(i)||s and by chain rule we have 

^WK*)II2) = IK*)II2 + ^IK*)II2- 

The second term in the above formula can be expressed as follows 

A A^v, A^v)h =2t(^v, A'v) H = —2t(Au, A s v) H + 2i(G 1 (i;), A s v) H 
at at 

+ 2t{G 2 (v + Z{6 t u)),A s v) H . 
Notice that for e > there are constants eg — cj.4 such that 
(J(rj,q),A s v) H < c 9 \\v\\ v \\r]\\ 1+s \\v\\ 1+s < c w (e)\\rj\\i +s \\v\\y + 

(J(ij,A v ),A s 

(J(r],Ar]),A s v) H < ciallA^Hill^lli+J-ylli+s < ci 4 (e)||»7||£||77||? +s + 



and 



J IMIs^ s — c i5 / IMIy^ s ) s < 1 



for the embedding constant C15 between 7i s and V. Then by the similar 
arguments as in P or [T3] we have the following dissipative property of 
THC. 

Lemma 4.4. For the random dynamical system if generated by hltfy . there 
exists a compact random set B = {5(a;)} CJGr2 which satisfies Definition 4-1- 

We define 



(28) B(u) = cp(l, 0-xlo, B(Q, R(0-iu))) CH S , < s < -. 

In particular, W is compactly embedded in H. 

Now we show that the dynamics of THC is described by a random attractor. 
First we recall the following basic concept; See [Zj. 
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Definition 4.5. Let tp be a random dynamical system. A random set 
A = {A(a;)} a ; g Q consisting of compact nonempty sets A(u) is called ran- 
dom global attractor if for any random bounded set D we have for the limit 
in probability 



for any t > and u G f2. 

The following theorem [2j gives a condition of the existence of random at- 
tractor. 

Theorem 4.6. Let if be a random dynamical system on the state space H 
which is a separable Banach space such that x — > (p(t,u,x) is continuous. 
Suppose that B is a set ensuring the dissipativity given in Definition (4-1)- 
In addition, B has a subexponential growth (see Definition (4-1)) and is 
regular (compact). Then the dynamical system ip has a random attractor. 

Then from the above analysis we have the following result to the random 
system <p generated by (jIBj). and, via the transformation (17), thus to the 
original stochastic THC model. 

Theorem 4.7. (Ttandom Attractor ) The THC model (l)-(3) has a ran- 
dom attractor in the phase space H = Lz(D) x L 2 {D) x L2(D). 

It is well known that, the global attractor for certain deterministic infinite- 
dimensional systems has finite dimension This result leads to the fact 
that asymptotic behavior of these systems can be described using finite- 
dimensional systems. And the similar theory has been developed for random 
dynamical systems; see, for instance, [51 ITU 123]. 

However, for the THC system we will apply another approach to prove the 
fact that the random attractor of (l)-(3) has only finitely many degrees 
of freedom. Namely we will use the concept of determining functionals in 
probability as in Chueshov et al. [$\. This concept was introduced by Foias 
and Prodi ^3] for deterministic systems. 

Definition 4.8. We call a set C — j — 1, ■ • • , iV} of linear continuous 
and linearly independent functionals on a space X continuously embedded 
in H (for instance X = 7i s orV ) asymptotically determining in probability 



(P) lim dist H (<f(t,cu,D(co)), A{6 t uj)) = 



and 



cp(t,u,A(u)) = A(O t u) 



if 




for two initial conditions v±, v?, G H implies 

(P) lim \\<p(t, u, v x ) - <p(t, oo, v 2 ) \\ H = 0. 
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We introduce a constant ec to describe the qualitative difference of the 
space H and X for some set of functionals 

(29) \\u\\ H < Crmaxl^u)! + £c|MU, C c > 0. 

We cite the following theorem from [Bj. 

Theorem 4.9. Lei £ = {lj : j = 1,...,N} be a set of linear continuous 
and linearly independent functionals on X . We assume that we have an 
absorbing and forward invariant set B in X such that for sup^^-, \\v\\ 2 x 
the expressiont — > sup„ gB ( etW ) \\v\\x is locally integrable and subexponentially 
growing. Suppose there exist a constant Ciq > and a measurable function 
I > such that for v\, v 2 G V we have for G(u, v) = Gi(v) + G 2 (v + Z(uj)) 

(-A(v x - v 2 ) + G(u, - G(cjj, v 2 ),v 1 - v 2 ) 

< -ciellfi - v 2 \\v + K^i^w) IK ~ v 2\\h- 

Assume that 




for some m > 0. Then C is a set of asymptotically determining functionals 
in probability for random dynamical system <p. 



Now we go back to the random THC system. Note that for e > there are 
constants cn(e), ■ ■ ■ , c 22 (e) such that 

\(J(ipi,qx) - J(i> 2 ,q2),qi -Q2)\ = \{J{$2-$iiqi),q2 - qi)\ 

<|||?i - h\\v + Ci7(£)\\qi\\%i {D )Mi - Q2\\l r 
Similarly 

s,) - j$ 2 , s*)s x - s 2 )\ < ||h - v 2 \\ 2 v + 019(6)11^11^11^ - s 2 \\l 2 , 

\(a T T y - a s S y ,qi- q 2 )\ < ^\\v x - v 2 \\\, + c w {e)\\q x - g 2 ||| 2 , 

and 

\(J($x,Ar}) - J(^ 2 , Ary),gi - q 2 )\ < c 21 \\vi - v 2 \\ 2 v \\r]\\ 2 w ,. 
Then we have 

(-A(vx - v 2 ) + G(w, vx) - G(u, v 2 ),vt - v 2 ) 
< -c 22 \\vi - v 2 \\\ + l{v 1 ,v 2 ,u)\\v 1 - v 2 \\ 2 H . 
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Here we can take k is large enough and e is small enough such that Ec 2 2 = 
1 - £ - c 21 E\\r]\\ 2 w ^ > 0. And 

l(vi,v 2 ,u) = ci 7 (£)||gi||^i (D) +ci 8 (e)||Ti||^i (D) + ci 9 (£)||5i||^ 1(D) +c 20 (e). 

Now we set X = 71 s , s E (0, |). In the above discussion we have shown 
that the set B, consisting of bounded sets, is forward invariant. Then we 
can apply Theorem 4.9 above to the random dynamical system (p generated 
by (|18|). and, via the transformation (17), get the following result to the 
original stochastic THC model. 

Theorem 4.10. (Finite Degrees of Freedom^ The THC model (l)-(3) 
has finitely many asymptotic degrees of freedom, in the sense of having a 
finite set of linearly independent continuous functionals which is asymptot- 
ically determining in probability, on 7i s with < s < 1/4. 

5. Conclusion 

We have investigated the dynamical behavior of a random thermohaline 
circulation model. We have shown that the random THC model is asymp- 
totically described by a random attractor ( Theorem 14. 7j) . And this system 
has finite degree of freedom in the sense of having a finite set of determining 
functionals in probability ( Theorem I4.10j) . 

Acknowledgement. We would like to thank Hongjun Gao and Bjorn 
Schmalfuss for helpful discussions. 
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